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00 ' Abstract 

o ■ 

■ Pade approximants are used to improve the convergence behavior of perturbative results 



in massless scalar and gauge field theories at finite temperature. 



iu: 

r-| , In recent years, computational methods have been developed to analytically tackle three-loop 
vacuum diagrams and higher-order contributions of diagrams with less loops in massless field theories 
at finite temperature [0-0- Consequently, the free energy density F at zero chemical potential could 
be computed analytically at the g^ level in both massless g'^(j)^ theory |^ (the pressure given there is 
the negative of the free energy density) and in massless gauge theories [0, In [|l-0, specializations 
to QED can be found, where the result was known before in partially numerical form 0. 

However, for interesting values of the coupling constant in non-Abelian gauge theories, the con- 
vergence behavior of the perturbative series is not convincing ^. In this brief report, we note that 
the use of Pade approximants drastically improves this behavior in both 0^ and gauge theories. For 
the use of Pade approximants and other resummation techniques in other contexts in perturbative 
field theory and statistical physics, see, e.g., [§, and references therein. 

Let us first review those features of the results of [0, ^, § which are essential for our analysis 
here. The perturbative series for the free energy density in both scalar and gauge theories has the 
structure (see the appendix for details) 

F = T^[C0 + C2g^ + C^g^ + (C4aln^ + C4b)g^ + (Csaln^ + C5fe)/ + O(^^ln^)] , (1) 

where T is the temperature and Cq, C2, C3, c^a, C5a are constants, while c^b and c^b have a logarithmic 
dependence on ln(/i/T), where ft is the renormalization scale in the modified minimal subtraction 
scheme (MS). In 0^ theory, = 0, while in gauge theories csa = 0. 
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As in 0, we use the renormalization group to make running, 

^-^-/5oln| + |ln(l-/5ogln|) , (2) 

where Pq and /?i are the one- and two-loop coefficients of the beta function Pg of g"^ (see the appendix 
for details on Pq and Pi), and gx is the coupling constant at temperature T. Then g in (|l|) is 
replaced by g{p). In this way we get an idea of the dependence of our result on the choice of 
renormalization scale. We could subsequently expand F in powers of gr to check that F becomes 
explicitly independent of ft through gj^. For this purpose, we would really only need the one-loop 
coefficient of Pg. The reason is that Pg contains only even powers of g, Pg = Pog'^ + Pig^ + ■ ■ ■, since 
we renormalize as at zero temperature. Therefore, from the viewpoint of the renormalization group, 
the g'^ and g^ terms in F are the ffist corrections to the g^ and g"^ terms, respectively. Numerically, 
the difference between using Pg to one or two loops is insignificant for the examples in non-Abelian 
gauge theories below, but keeping the two-loop correction turns out to improve the behavior of the 
resummation in 0^ theory. The reason why it suffices to use the two-loop beta function in the 
examples in this work is that bad behavior of both the perturbative result and Fade approximants 
sets in for values of g^ where the two-loop beta function is still a good approximation. 

Now we use Fade approximants to reexpress F. For this purpose we pretend that C4 = C4a In (7+04;, 
and C5 = csa hig + c^b are constants in a Taylor series F = T'^[cq + C2g^ + c^g^ + c^g^ + c^g^ + ■■■]. 
C4 and C5 have different values for each choice of /2 both through their direct dependence on In (/i/T) 
and through the running of g{ji) in \ag. Using the approximants [1/2], [2/2] and [2/3] to rewrite F 
through orders g^ , g^ and g^ (there is no approximant [1/1], since F contains no term linear in g) 
gives 

TP, , — 7^4 cgc2-cgc3g /q\ 

7-1 _ rpiCoc\-CQC2Cig+{c\+COci-COC2C4)g'^ ( ,\ 

^[2.2] - ^ ci-C2Czg+{ci-clc^)9^ ' W 

p _ rp4: Co{c^+Cocl-C^C2C4)+CQ(-clc3-CoC3C4+CoC2C5)g+(c^+2coC2cl-^ /g\ 



Define 



a{T) = g^iAn) (6) 

and let us look at some specific examples. Our ffist case is the small- coupling QCD example from 
[01 with dA = 8, Ca = S, Uf = 6, dp = 18, Sp = 3, S2F = 4, and a(T) = 0.001. As argued in 
[^ and as can be seen in Fig. la, the perturbative series for F through g^ is well behaved in this 
case, with respect to both convergence for a given renormalization scale ft and to the growing ft 
independence of F towards higher orders. The Fade approximants -^[1,2] and -^[2,2] are close to the g^ 
and g'^ results within the expected accuracy (given by the magnitude of the g^ and g^ corrections, 
respectively). However, F[2,3] has a pole, as seen in Fig. lb. This pole comes about through a zero of 
the denominator in (0), which, in turn, due to the smallness of a{T), is caused by a zero of the ffist 
term in the denominator of (^, c|-|-Coc| — C0C2C4. We know that the full result for F is independent of 
ft and that consequently this pole is an artifact of the resummation scheme. We therefore determine 
its position and residue, explicitly remove it and call the resulting function -^[2,3]- The curve in Fig. lb 
for -F[2,3] is virtually identical to the g^ result in the perturbative series in Fig. la. 



2 



Now let us turn to cases where the pure perturbative series needs improvement. Fig. 2a represents 
the perturbative series for the pure SU(2) example from |^ with rf^ = 3, Ca = 2, nj = dp = Sp = 
S2F = 0, and a(T) = 0.03, while Fig. 2b shows the Fade approximants. Again, we have removed 
the pole from -F[2,3] to define -^[2,3] and show both functions. Clearly, the convergence behavior 
of the series -^[1,2], ^[2,2], -^[2,3] is drastically improved compared to the purely perturbative series, 
particularly around natural choices of fl, such as /2 = T or /i = gT, which up to order are the only 
mass scales in finite-temperature non-Abelian gauge theories. Note also the relative independence 
from the renormalization scale. 

Now turn to the other QCD example in [^], namely, a(T) = 0.1 with the other parameters being 
the same as in our first case. The result is plotted in Fig. 3. Again, there is much improvement 
compared to the pure perturbative series around fi = T ^ gT, where higher approximants give 
subsequently smaller corrections to their predecessors. 

As our final example in non-Abelian gauge theories, consider three-flavor QCD, i.e., = 8, 
Ca = 3, Uf = 3, dp = 9, Sp = 3/2, 5*2^ = 2 with a{27iT) = 1/3 [note that we have to replace 
T — > 2ttT in (|]) and (H) accordingly]. Up to the fact that we neglect the strange-quark mass 
and that we have set all chemical potentials to zero, this is close to the case of the quark-gluon 
plasma to be produced at the BNL Relativistic Heavy- Ion Collider (RHIC). The result is plotted in 
Fig. 4. There seems to be no useful improvement over the perturbative series, although the range of 
manifestly bad behavior is shifted towards smaller values of /2. 

In Fig. 5, we present an example in scalar theory, namely with a{T) = 0.75. Note how, at least 
for not too big /i, the Fade approximants fluctuate much less in subsequent orders than their purely 
perturbative counterparts. The fact that we can go to larger couplings in scalar theory than in 
non-Abelian gauge theory is easily explained. For example, for the case of no fermions, the effective 
expansion parameter in F is seen to be Ca«(T). That is, a larger number of degrees of freedom 
leads to stronger corrections to the ideal gas result (unless we try to make fermionic and bosonic 
contributions cancel). 

Let us make two final remarks, (i) The use of the approximants [2/1] and [3/2] instead of [1/2] and 
[2/3] gives results very similar to those presented here, while approximants [m/n] with \m — n\ > 1 
give typically less improvement, (ii) Starting at g^ order, another physical scale g^T enters the 
calculation of F in non-Abelian gauge theories [^. Therefore, it would be interesting to see how 
inclusion of the g^ term changes our results. Unfortunately, computation of this term is difficult and 
requires a combination of perturbative and nonperturbative techniques |T0 . 
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Appendix 

Here we give the results of |^ and ||, ^. 
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With the Euchdean Lagrange density 



(7) 



the free energy density in the MS scheme is 
10 8 V47ry V6 V47ry 
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where we have translated the MS resuh of ||^ into MS using /x^ = e^^ ^ j (47r) and where C is Riemann's 
zeta function and 7e is the Euler-Mascheroni constant. The one- and two- loop coefficients in /3g are 



/3o 



(47r)^ 



17 



3(47r) 



(9) 



In gauge theory with fermions with a single, simple Lie group consider the Euclidean Lagrange 
density 

^ = i^l, (d, - igA^T^) ^ + ^ [d^Al - d^Al + gt^'AlA'^^ , (10) 

where the T°- are the generators of the group in the fermion representation. Let (Ia and Ca be the 
dimension and quadratic Casimir invariant of the adjoint representation, with 



ahc rdbc 



dA, nr 



:ii^ 



Let dp be the dimension of the total fermion representation {e.g., 18 for six-flavor QCD), and deflne 
Sp and 5*2^ in terms of the generators T*^ for the total fermion representation as 



Sf = — tr(T2) , = — tr[(T^ 

cLa cLa 



'2\2i 



(12) 



where = T^-T"". For SU(A^) with np fermions in the fundamental representation, the standard 
normalization of the coupling gives 



- 1 , Ca = N , dF = Nnp 



1 N"^ -I 

Sf = -np , S2F = — — np . 



(13) 



The the free energy density is given by 
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The one- and two-loop coefficients in j3g are 
1/22 
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Figure 1: Fig. la shows the perturbative series for the free energy density in units of the ideal gas 
result F{g = 0) for six-flavor QCD with a{T) = 0.001 for a range of choices of renormalization scale 
p,. The short-dashed, medium-dashed, long-dashed and sohd lines are the results for F including 
terms through orders g^, g^, g'^ and g^, respectively. Fig. lb shows Fade approximants instead: The 
g'^ result has been dropped, while the result through g"^ from Fig. la has been replaced by -^[1/2], the 
result through g^ by -Fp/2] and the result through g^ by -^[2/3] (solid line with pole) and -Fp/3] (solid 
line without pole). 
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Figure 2: The same as Fig. 1, but for pure SU(2) theory with a{T) = 0.03. 
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Figure 3: The same as Fig. 1, but with a{T) — 0.1. 



9 



m 



m 



1.2 



1 . 1 



1 - 



0.9 



0.8 



1.2 



1 . 1 



0.9 



0.8 



■ \ 

K 






/ / 
/ / 
/ / 
// 




2 


4 6 8 

logio/i/T 
Fig. 4a 












/ ^ 

/ 

/ 
.1 




2 


4 6 8 

logio/^/T 

Fig. 4b 



Figure 4: The same as Fig. 1, but for only three fermion flavors and with Q;(27rT) = 1/3. 
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Figure 5: The same as Fig. 1, but in scalar theory with q;(T') = 0.75. 
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